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The ground and excited states of *He were investigated with a method of antisymmetrized molec- 
ular dynamics (AMD). We adopted effective nuclear interactions which systematically reproduce the 
binding energies of *He, ^He and *He. The ground state of ^He has both the j-j coupling feature(p3/2 
closure) and the L-S coupling feature(''He-|-2n + 2n) with a slight tail of dineutron at the long dis- 
tance region. The theoretical results give an indication of the 0^ state with dineutron gas-like 
structure. The dineutron structure, '^}ie+2n+2n, of this state is similar to the 3a-cluster struc- 
ture of the ^^C(Oj) state which has been interpreted as an a condensate state. Since the *H[e(Oj) 
state has a significant overlap with the dineutron condensate wave function where two dineutrons 
are moving in S wave around the a core with a dilute density, we suggest that this theoretically 
■ predicted state is a candidate of the dineutron condensate state. 
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I. INTRODUCTION 

In the recent progress of unstable nuclear physics, various kinds of exotic structure have been discovered. Many of 
these phenomena in light nuclear region are often related to cluster physics. From the viewpoints of nuclear cluster, 
there are many theoretical works on halo structure in neutron- rich nuclei and molecular structure in Be isotopes. 
Recently, Tohsaki et al. proposed a new type of cluster structure in the second 0+ state ^^C, where 3 a clusters are 
weakly interacting[l|. This is a dilute gas state of a particles which behave as bosonic particles in the dilute density. 
Q , This phenomena is associated with Bose-Einstein Condensation(BEC) and is called "alpha condensation" . The alpha 
^ ' condensation was originally suggested in dilute nuclear matter by Ropke et al.Q. The 0^ of ^^C is regarded as an 
(~| [ example, where the alpha condensation is realized in a finite nuclear system. Then, it is challenging to search for such 
cluster-gas states in other nuclei. In analogy to the alpha condensation, dineutron condensation in neutron matter 
^ ' is a recent key issue in physics of unstable nuclei. Matsuo suggested that the dineutron correlation can be enhanced 
^ in dilute neutron matter [3|. In real systems, one should focus on dineutron correlation in finite nuclei such as halo 
nuclei and extremely neutron-rich nuclei, or that in neutron skin at a surface region of neutron-rich nuclei. In fact, 
CS| , the dineutron correlation in two-neutron halo nuclei like ^He and "'^^Li attracts great interests in these days. In case 

■ of ^He, where the *He is the good core, the dineutron correlation of the valence neutrons has been demonstrated in 
three-body model calculations (for example, [4, i^, [1,01 and references therein). 

[ — ■ Now, let us consider structure of ^He from a point of view associated with the dineutron condensation. Firstly, 
more than one dineutrons are required to construct a dineutron condensate state. In *He, two pairs of neutrons are 
C possible from four valence neutrons around the ^He core. In second, ®He system may have some correspondence with 

■ the ^^C system, because both of them have the same neutron number, = 6. In analogy to ^^C, the ground state 
J> . of ®IIe may have a feature of the neutron P3/2 closure or the SU(3)-limit p-shell configuration. Instead of the ground 

state, one can speculate the dineutron gas-like state with developed '^Re+2n+2n structure in excited states. 
^ There are many theoretical works on He isotopes. Application of ab initio calculations such as GFMC and NCSM 

d ' with realistic nuclear forces have now reached to the mass A ~ 10 region including ^He and ®He[l, [^[l^. Systematic 
studies of He isotopes have been performed also by model calculations with effective interactions such as cluster models 



as well as GSM[li|, [l^jliSl and mean field a ppr oaches 14 1. Three-body model with an assumption of the ''He core has 
been often adopted to study ^He [1, SIEHllEIlll and it has been a ppli ed to heavier He isotopes JTj]. ®He and ^"He 
have been also studied by such models as ^He-|-Xn models [3 fisl l20l [HI and by extended models [13, [2^ which 
have less assumption of the core. With Fermionic molecular dynamics, the study of He isotopes has been performed 
based on a realistic nuclear force [24] ■ However, many of these studies are concentrated on the ground states except 
for three-body models, GSM and GFMC. 

After the experimental indication of neutron skin structure in ^He[25j, many experimental works on ®He have 
been recently performed to reveal the detailed properties of the ground state. The core excitation ^He(2+) in the 
ground state, which has been experimentally suggested (26|. indicates that *He is different from a simple three-body 
state of ^He(0"'")-|-2n. Recent experiments using ®He beams suggested the significant component of the (^3/2)^(^1/2)^ 
configuration [23, [ll]. They may support dineutron correlation in the ^He ground state rather than the pure (^3/2) 
closure of neutrons. On the other hand, a measurement of spectroscopic factor of ''He(3/2~)[2^ in ^He suggested the 
pure sub-shell closed structure contradictory to the other experimental results. Thus, the neutron structure of the 
®He ground state is controversial. Concerning excited states, although some levels are known to exist in the energy 
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Ex = 3 ^ 8 MeV region, the experimental information is very poor for these states except for the 2^^ state [3C|. 

In this paper, we investigated structure of ^He. In particular, we focused on 0+ states and discuss their dineutron 
component, because one of our major aims is to search for the dineutron gas-like state. We applied a method of 
antisymmetrized molecular dvnamicsfAMD)[3ll. [3^ . [ssj . which has been already proved to be useful in describing 
cluster structure in light nuclei. AMD has been applied to various light unstable nuclei such as He, Li, Be isotopes 
as well as stable nuclei. It has been applied also for study of cluster gas-like states in ^^C and ^^C(^^B)[33. l35l |. 
In the present work, we adopted a AMD-|-generator coordinate method(GCM). Namely, we superposed a number of 
AMD wave functions, which were obtained by energy variation with constraints, to take various configurations into 
account. We comment that the theoretical method AMD-I-GCM of the present calculation is similar to those of the 
AMD-t-GCM and AMD-I-SSS works on He isotopes by Itagaki and his collaborators [13, [13| in a sense that multi 
configurations of AMD wave functions are superposed. In [2Ct [23j , ''He-I-Xn and t + t + Xn configurations were a 
priori assumed. Another claim is that they used an effective interaction which makes a bound ^n. In the present 
work, we have no assumption of the cluster core and chose effective interactions by taking care of subsystem energies 
such as a-n and ^He as well as nucleon-nucleon scattering. We used some sets of interaction parameters and showed 
the calculated results of the ground and excited states of He isotopes. By assuming (Os)^ configuration as the interior 
structure of a dineutron, we analyzed dineutron structure of *He and compared it with the a-cluster structure of ^^C. 

The paper is organized as follows. In the next section, we briefly explain the theoretical method of the present 
work. Results are given in IIIIl and dineutron structure is discussed in IIVI Finally, we give a summary in |Vl 



II. FORMULATION 



In this section, we briefly explain the formulation of AMD-I-GCM in t he p resent calculation. The detailed formu- 
lation of the AMD method for nuclear structure study is described in [s^, Hj] . There are various versions of practical 
methods of the AMD framework. In the present work, we performed superposition of a number of AMD wave func- 
tions obtained by energy variation with constraints based on the concept of GCM. The procedure of the variation, 
spin and parity projection and superposition is similar to those of AMD-I-GCM calculations in [13, [13, [37i], though 
the details of model wave functions and effective interactions are different from each other. 

An AMD wave function is a Slater determinant of Gaussian wave packets; 

<^AMd{'Z) = -^A{ipi,ip2,---,'PA}, (1) 

where the ith single-particle wave function is written by a product of spatial((/)), intrinsic spin(x) and isospin(r) wave 
functions as, 

= (j^yi^Xin, (2) 
0x.(r,) = (^|yexp{-Kr.-^r}, (3) 

X^ = {\+idx^ + {\-^i)Xl- (4) 

(/)Xi and Xi are spatial and spin functions, and Tj is isospin function which is fixed to be up(proton) or down (neutron). 
The width parameter u is chosen to be the optimum value for each system. Accordingly, an AMD wave function is 
expressed by a set of variational parameters, Z = {Xi, X2, • • • , X^, ^1, ^2, • • • , Ca}- 

The energy variation was performed for the parity-projected AMD wave function <I>^j^j-,(Z) under constraints. In 
order to obtain basis wave functions, we adopted the total oscillator quanta and deformation as the constraints. 
Hereafter, we note the expectation value of an operator O with respect to a normalized parity-projected AMD wave 
function as (O). Expectation values {N^°) of the total oscillator quanta is given by the creation and annihilation 
operators of harmonic oscillator in the same way as [s^l- In the AMD-I-GCM calculations with the /3-constraint (for 
example [131), the deformation is usually constrained by using the rotational invariant value D = Tr{QQ)/Tr'^{Q), 
where the matrix Q is calculated by quadrupole operators as Qa-p = ^iPi) (o" = x,y,z and p = x,y, z) [38|]. Here 
D is approximately related to the quadrupole deformation parameter /3 as l?(/3) = (5/3^/27r -I- l)/3. In the present 
work, we used the modified quadrupole matrix Q'^^ = Qa-p — AS^p {A is the mass number) instead of the original 
Qap and imposed the constraint on the D' = Tr{Q'Q')/Tr^{Q'). This is useful for He isotopes to obtain basis wave 
functions with various configurations on mesh points of the two-dimensional parameters, (3 and {N^°). The energy 
variation with the constraint values A'const and /3const was performed with respect to the parity-projected AMD wave 
function by minimizing the energy defined as, 

E={H)+ V^{N,onst - (/^'°))' + V^^(i^(/3const) " D'f. (5) 
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Here the artificial potentials are introduced to satisfy the condition of the constraints. With a given set of constraint 
values (iVconst) /Sconst) the optimum wave function $amd (-^const j /Sconst ) was obtained. Finally, we superposed the 
spin-parity eigen states projected from the obtained wave functions, 

|«He(J±))= C;^^(^const,/3const)|PA^/;f$AMD(^const,/3co„st)), (6) 

JVco„st,/9con=t 

where the coefficients c^^ A'const , /Sconst) were determined by diagonalizing the Hamiltonian and Norm matrices. In 
the present calculations, we took only M = K = states. 

III. RESULTS 
A. Calculations 

^He, *He and ^"He were calculated by the AMD+GCM method. The strengths, and V^^, for the constraint 
potentials in cql5] are chosen to be 30 MeV and 2000 McV, respectively. We chose the width parameter v to opti- 
mize the energy for the ^'('j\7/c)=(oo)^AMD(-^const = -^min + 2), which gives the minimum energy among the states 
^('MA')=(oo)^AMD(^const) in most cascs. Here, A^n^in is the minimum value of the harmonic-oscillator quanta, N^^i^ =2, 
4, and 6 for ^He, *He, and ^°He, respectively. A common value for each He isotope are used in the calculation 
with each interaction. The adopted v values are listed in table HI We adopted the constraint values of the mesh 
points on the A^const-/3const plane as A^^oLt — ^min + A(*)(A(*) =0,1,2,3,4,6,8,10 for positive parity states and 

A('^ = 1,2,3,4,6,8,10 for negative parity states) and /3Hnst=0; 0-2, 0.4, 0.6, • • •, 1.6. Then, the total number of the 
basis wave functions are 72(63) for positive(negative)-parity states. On the A"const-/?cont plane, we first obtained the 
wave function <i>AMD(-^const, /3const) at A^const = A^min + 2 and /3const=0, 0.2, 0.4, 0.6, • • •, 1.6. Then we searched 
for $AMD(-^const + l,/?const) (or ^AMol^const " l,/3const)) Starting from the $AMD(^const,/5const) by increasing(or 
decreasing) A^const one by one. 

Some of the basis wave functions with the constraints have the breaking of the *He-core. Such the basis wave 
functions with the ^He-core breaking have high energies in general, and therefore, they practically give only small 
contribution to the low-lying states of ^He, ^He and ^°He isotopes. It means that the ''He cluster is a rather good 
core in ^He, ^He and '°He isotopes, while the motion of valence neutrons is relatively important. 

TABLE I: Parameter sets of the effective interaction and the values of width parameter v adopted in the present work. The 
theoretical values of scattering length as{at) for singlet (triplet) even channel, neutron separation energy of ^He (5'n(^He) = 
£;(''He)-£'(''He-n)), 2a threshold energy of *Be, two-neutron separation energies of ''He and *He {S2nfVie)= E{'^Re)-E{^ne) 
and 5'2n(**He)= £;(''He)-£;(*He)) are also listed. 



Parameter set 




v58 


v56 


m62 


m56 


Central force 




Volkov No. 2 


Volkov No.2 


MVl case(3) 


MVl case(3) 


Wigner 


w 


0.42 


0.44 


0.38 


0.44 


Bartlett 


b 





0.15 





0.15 


Heisenberg 


h 





0.15 





0.15 


Majorana 


m 


0.58 


0.56 


0.62 


0.56 


i^(^He) (fm-2) 




0.265 


0.265 


0.210 


0.210 


i^i^Ue) (fm-2) 




0.245 


0.245 


0.210 


0.210 


:/(»He) (fm-2) 




0.240 


0.240 


0.185 


0.185 


vC°Re) (fm-2) 




0.185 


0.175 


0.165 


0.165 




exp. 


v58 


v56 


m62 


m56 


at (fm) 


5.42 (p-n) 


9.7 


5.4 


6.4 


4.2 


as (fm) 


— 16.5 (n-n) 


9.7 


-23.9 


6.4 


>100 


SnCRe) (MeV) 


-0.9 


-0.7 


-0.7 


-1.0 


-0.4 


2_E(^He)-£(''He-''He) (MeV) 


-0.1 


0.6 


1.4 


-1.3 


-0.6 


S'2n(®He) (MeV) 


1.0 


1.3 


-0.2 


2.1 


1.1 


S2nCHe) (MeV) 


2.1 


3.0 


3.2 


1.2 


2.0 
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B. Interactions 

We used effective nuclear interaction consisting of the central force, the spin-orbit force and Coulomb force. As for 
the central force, we adopted the Volkov force|39!] used in the work on He isotopes with AMD-l-GCM(^He+Xn')[20j. 
and also the MVl force [4^ used in the AMD calculations of ^^C [13, SI]. We used the spin-orbit force of the G3RS 
force fill as done in [20l. |42{ . We fixed the strengths of the spin-orbit term as uis = 2000 MeV, which is the same value 
as in [20]. By taking care of energies of subsystems, we tuned the interaction parameters, w, b, h, m, for Wigner, 
Bartlett, Heisenberg and Majorana exchange terms in the the central force(Volkov or MVl), respectively. ^He, ^He 
and ^°He were calculated with AMD-I-GCM by using totally 4 cases of central force. The parametrization for the 
central force is summarized in table HI In order to demonstrate characteristics of the effective interactions, we also 
show the relative energies of subsystems and the nucleon-nucleon scattering lengths with these 4 types of interaction. 
We estimate the energies of the ''He, ''He-n state with J'^ — 3/2^, and ■*He-'*He state with =0+, by assuming the 
(Os)^ state of *He and performing cluster-GCM calculations within the a-n and a-a cluster models for simplicity. 

The first case of interaction is Volkov No. 2 force [39| with interaction parameters m = 0.58, b — h = 0. This is 
the same effective interaction as that used in the AMD-l-GCM(^He-l-Xn) by Itagaki et al. [13], which succeeded to 
systematically reproduce the binding energies of He isotopes. We note this interaction 'v58' in this paper. In spite of 
good agreement of the binding energies of He isotopes, the v58 force has a fault that 2 neutrons are bound in a free 
space. It is well known that the Volkov force with b — h = Q has too strong neutron-neutron attraction, because such 
the parametrization with no Bartlett term nor Heisenberg term gives the same interaction in the singlet-even channel 
as that in the triplet-even channel. In reality, the singlet-even channel has weaker attraction, and two neutrons are 
unbound. In order to describe dineutron correlation in neutron-rich nuclei it might be crucial to reproduce such the 
feature of two-nucleon system, though it does not matter in case of spin-isospin saturated systems like Z = N nuclei. 

In the second case of interaction, we used Volkov No. 2 force with modified interaction parameters as m = 0.56, 
b = h — 0.15. This interaction (noted as 'v56') describes well the experimental 5'-wave scattering lengths of the n-n 
and p-n channels, and the unbound feature of 2-neutron system. The Majorana parameter m — 0.56 was determined 
by adjusting the binding energy of ^He to the experimental data. However, this interaction fails to reproduce 2n 
separation energies of ^He and ^He, and it also gives too strong attraction in '*He-'*He system. 

The third interaction('m62') and the forth one('m56') listed in table U are based on the MVl force[40i|. The 
parametrization of the m62 interaction is m = 0.62 and b = h = 0, which is the same as used in the AMD calculations 
of ^^C [sl-li^. In case of the m62 interaction, two neutrons are bound in a free space as well as the Volkov force with 
b = h = like the v58 interaction. In the 'm56' interaction, we used the modified Bartlett and Heisenberg terms, 
6 = /i = 0.15, and the Majorana term m = 0.56 which was adjusted to reproduce the binding energy of ®He. With 
the TO = 0.56 interaction, two neutrons are almost unbound in a free space, and other energies of subsystems are 
reasonably reproduced. 

C. Ground states of He isotopes 

We show the calculated results of the ground states of He isotopes. The energies of He isotopes arc shown in Fig. [T] 
The v58 and m56 interactions systematically reproduce the energies of "^He, ^He and ^He, though they overestimate 
the ^'^He energy. On the other hand, the v56 and m62 interactions are poor in reproduction of the ^He energy, and 
therefore, they fail to reproduce two-neutron separation energies of ^He and *He as shown in table [H Hereafter, 
we discuss the results obtained with the v58 and m56 interactions. We stress again that the v58 interaction well 
describes the energies of subsystems except for the fault of the too strong neutron-neutron interaction, while the m56 
interaction reasonably reproduces the global features of the subsystem energies. 

The calculated root-mean-square radii of proton, neutron and matter density are given in table [H] with the ex- 
perimental data. The theoretical results of other calculations are also listed. Experimentally, extremely large radii 
of ^He and ®He have been reported by the reaction cross sections [1^, 13] ■ It has been suggested that the large 
radii originate in the remarkable enhancement of neutron radii due to the neutron-halo and neutron-skin structures 
in ''He and ^He, respectively. The empirical neutron radii are well described by the present calculations with the 
m56 interaction. On the other hand, the neutron radii calculated with the v58 interaction are slightly smaller than 
the empirical ones as well as the former AMD-l-GCM(^He-l-Xn) calculations with the same v58 interaction[23]. The 
proton radii calculated with the m56 interaction are consistent with the observed data except for that of ^He. Figure 
[5] shows the proton density and neutron density. In ^He, the neutron density has a long tail at a large distance 
region. This is the neutron halo structure and is similar to the neutron density obtained by other calculations such 
as SVM[19|. In *He, the neutron and proton density shows the neutron skin structure at the surface, which well 
corresponds to the discussion in [l9l, [23. Thus, the present calculations with the m56 interaction systematically 
describe the ground-state properties of ^He and ®He such as energies and radii. 
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FIG. 1: The calculated energies of He isotopes with the v58, v56, m62 and m56 interactions(see text). The experimental data 
are also given. 



Let US discuss the effect of the spin-orbit force, which may induce the j-j coupling feature of neutrons. The 
expectation values of the spin-orbit force (Vis) and those of the squared total intrinsic spin of neutrons (S^) are listed 
in table Hill From the values of (S*^), the S = 1 component in the ^He(0]^) state is estimated to be 0.13 and 0.07 
in the m56 and v58 results, respectively. It means that the (7*3/2)^ configuration is contained due to the spin-orbit 
force. However, the 5 = component is still significant because of L-S coupling feature of spin-zero 2n correlation. 
We note that the fraction 0.87 in the m56 results for the 5 = component in ^He is in good agreement with three- 
body model calculations 0, [H, [3, HE] • Compared with the results of ^Hc, where the L-S coupling configuration is 
significant as well as the j-j coupling configuration, the j-j coupling feature increases in the *IIe(0]^) state because 
of the (^3/2)^ closure. As a result, the spin-orbit force gives much larger attraction in ^He by factor 3^4 than in 
^He. It is interesting that the the value (Sf-^) = 0.86(0.72) of the ^He(0^') in the m56(v58) results is different from the 
value (Sf-i) = 1.33 for the pure (P3/2)'* closed state. This deviation is because the L-S coupling configuration is still 
contained in ^He due to the spin-zero 2n correlation of neutron pairs. The detailed dineutron structure of ^He and 
^He will be discussed later. 

TABLE II: Root-mean-square radii (fm) of point-proton, point-neutron and point-matter density of the ground states of He 
isotopes. The experimental value(a) is deduced from the charge radius 46], and em piri cal values(b) are taken from [25l. |4^. 
Theoretical values of other calculations, NCSMlOl, SVM;U] AMD-(-GCM(''He-fXn)[23|, RMFl4i| are also given. 







exp. 


AMD-v 


58 AMD-m56 SVM[19] RMF[14] AMD(^He-HXn) 


[20] NCSM[10] 


*He 




1.455(1) 


1.46 


1.64 






1.45 




rn 




1.46 


1.64 






1.45 








1.46 


1.64 




1.76 




«He 


rp 


1.912(18)'"^ 


1.83 


1.90 


1.80 




1.89 




r„ 


2.59 - 2.&l'-^'> 


2.40 


2.49 


2.67 




2.67 






2.33 - 2.48*'') 


2.23 


2.31 


2.46 


2.43 2.32 




«He 


rp 


1.76 - 2.15*''' 


1.76 


1.96 


1.71 




1.88 




r„ 


2.64 - 2.m'-^'> 


2.37 


2.63 


2.53 




2.8 






2.49 - 2.52^'') 


2.24 


2.48 


2.40 


2.55 2.31 




i»He 


rp 




2.04 


2.13 










rn 




2.88 


2.97 










rra 




2.73 


2.82 




3.17 
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FIG. 2: Point-proton and point-neutron density in the ground states of He isotopes. The calculated results are those with the 
m56 and v58 interactions. 

D. Excited states of *He 

The calculated energy levels of *He are illustlated in Fig. [31 and the properties of the excited states are shown 
in table IIIII In both of the in56 and v58 results, the 2^ state is the lowest excited state and the state appears 
just above the 2f state. The lj~ and 3j~ states are obtained in a higher energy region. In addition, in the present 
calculations with the m56 interaction, the if , Oj~ and 2j~ states are obtained in almost the same energy region as 
the 1^ and 3j~ states. The present AMD framework is regarded as a kind of bound state approximation because of 
the restricted model space, and therefore, coupling with continuum states is not taken into account. In such a case, 
only resonance states remain in low-energy region while continuum states rise to a high excitation energy region in 
principle. However, in order to check the stability of the resonances against neutron decays, their properties should be 
carefully examined. In the present m56 results, the negative-parity states contain large component of ®Hc+n + rt-like 
configurations with the valence neutron far from the core. Since they have extremely large neutron radii and show 
somehow escaping behavior of neutrons, further investigation is required for these negative-parity states. In particular, 
the 2^ and 0^ states can couple with (Os)^(Op)'^(ls)^ neutron configuration which has a valence lsi/2 neutron 
with no centrifugal barrier. 

Compared with the experimental data, the theoretical values of the 2^ excitation energy are higher than the 
experimental one. However, it is important that the level structure for the excited states, , Oj, 1^ and 3^ , is not 
sensitive to the adopted interaction though the relative position to the ground energy depends on the interaction. 



7 



> 



15 
10 
5 

-5 



''He+4n.. 
f'He+2ri"" 



-2+ 



exp 



3" - 
2" I 
1+= 



•'2+ 



AMD 
m56 



-3" 

"r 

10+ 

2+ 



AMD 
v58 



FIG. 3: Energ 
taken from |47 



levels of He. The calculated results are those with the m56 and v58 interactions. The experimental data are 



TABLE HI: Excitation energies, Root-mean-square radii of point-proton, point-neutron and point-matter density, the expec- 
tation values of squared total intrinsic spin of neutrons (5^), and those of the spin-orbit force {Vis)- 







exp. 


AMD-V58 












AMD-m56 












nucleus 


Ex 


Ex 








(Si) 




Ex 


Tp 


r„ 




{Si) 


{Vi.) 






(MeV) 


(MeV) 


(fm) 


(fm) 


(fm) 




(MeV) 


(MeV) 


(fm) 


(fm) 


(fm) 




(MeV) 


«He 




1.797 


3.2 


1.82 


2.42 


2.23 


0.19 


-2.3 


2.6 


1.87 


2.46 


2.28 


0.27 


-2.3 


''He 


0^ 





0.0 


1.83 


2.40 


2.23 


0.16 


-2.6 


0.0 


1.90 


2.49 


2.31 


0.26 


-2.3 


«He 


or 
















10.8 


2.13 


3.63 


3.32 


2.05 


-5.9 


^He 


2r 
















10.8 


2.07 


3.41 


3.13 


2.00 


-6.2 


^He 


1^ 
















9.0 


1.94 


2.81 


2.62 


2.03 


-2.5 


^He 


3r 


7.16 


13.5 


1.90 


2.89 


2.68 


0.64 


-6.7 


11.5 


2.09 


3.31 


3.05 


1.02 


-5.3 


^He 


ir 


4.36 


12.1 


1.95 


3.05 


2.82 


0.81 


-7.9 


9.8 


2.13 


3.52 


3.23 


1.24 


-5.8 


«He 






10.3 


1.97 


2.94 


2.73 


0.67 


-4.7 


8.5 


2.11 


3.12 


2.90 


0.99 


-1.0 


^He 




3.1 


9.3 


1.76 


2.48 


2.32 


0.39 
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2.97 
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The 0^ state is theoretically suggested to appear just above the 2^ state. What is striking is that the state has 
a remarkably large neutron radius compared with the ground state because of developed ^Hc + 2n + 2n structure. 
In the obtained wave function of the 0^ state, which is given by a superposition of the basis AMD wave functions, 
the amplitude is found to be widely distributed into the basis wave functions with various spatial configuration of 
^He -f 2n -I- 2n. This indicates a gas-like feature that the dineutrons are rather freely moving around the ^He core. 
Therefore, we consider that the state is the candidate of the cluster gas-like state with two dineutrons around the 
a core. The detailed discussion of the dineutron-like structure is given later. In the experimental energy spectra, some 
excited states were observed above the 2^ state. Spins and parities of these states are not definitely assigned yet. In 
the present calculations, the predicted state has the strong monopole neutron transition from the ground states 
as the matrix element M„(05^ Oj) = 13.5(13.9) fm^ in the m56(v58) results. This neutron matrix element is much 
larger than the observed proton matrix element M„(0^ — > Oj) = 5.4 fm^ of ^^C by more than factor 2. Therefore, 
we consider that the ^He(Oj) might be excited by inelastic scattering on nuclear target. 

The excited states of ^He have been theoretically predicted by a few other calculations such as CSM and GFMC. 
The CSM gives better agreement of the 2^ excitation energy with the experimental data[l3]. We also comment 
that the GFMC calculation with AV18/IL2, which is an ah initio calculation with the realistic 2-body force and the 
empirical 3-body force, gives similar level structure to the present m56 results. Namely, the GFMC with AV18/IL2 
gives the 2+ state at — 4.72 MeV and the Ij^, and 2^ states in the > 5 MeV region. 
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IV. DINEUTRON STRUCTURE 
A. What is dineutron(^n) cluster ? 

There is no bound state in an isolate nn system. However, it has been emphasized in many theoretical works 
that the spatial neutron-neutron correlation plays an important in the binding mechanism of the Borromean systems 
with two-neutron halo such as ^He and ^^Li (for example, [1, SSSl] and references therein). The neutron-neutron 
correlation is characterized by a spin-zero nn pair with spatial correlation in S wave. In the correlation density of 
two-neutron halo nuclei, a peak of the probability appears at the region with a small n-n distance(i?(nn)) and a large 
n-core distance in general. This corresponds to the dineutron correlation. In an extended meaning, it is regarded as 
a "dineutron cluster" which can virtually exist in loosely bound neutron-rich nuclei. 

As mentioned above, the characteristics of the dineutron are the zero spin and the spatial correlation. In the 
correlation density for ^He, ^^Li and ^^Be given by three-body calculations 0, 0) IH) 113] > the peak for the dineutron 
correlation are seen typically around the R{nn) ~= 2 fm with a ridge in the R{nn) = 2 ~ 3 f m region. It is important 
that this n-n distance at the peak nearly depends on the system size among these three systems, ^He, ^^Li and ^''Be. 
From this most probable n-n distance, the typical size of the spatial correlation of the nn pair can be estimated to 
be about 2 fm. Then, we here approximately describe the dineutron cluster, ^n, by a spin-zero neutron pair written 
by the simple harmonic-oscillator (Os)^ state with the size parameter b in order to investigate dineutron structure in 
*He. Then, the ^n-cluster wave function cj) "(S) which is localized at the position S is expressed as, 

</>^"(S)=^{<^«^(ri)xT0"s^(r2)Xi}, (7) 
'/'^^(r.) = ^exp{-5l,(r,-S)2}. (8) 

In this definition, the relative motion between two neutrons in the ^n cluster is given by a Gaussian, 

</.'-(ri - r^) = — exp{--L(ri - r^)^}, (9) 

with the size br = \f2b^ which should be the typical nn distance br — 2 ^ Z fm. With this approximation of the 
cluster, major component of the dineutron correlation might be taken into account, though the tail part at the 
large correlation length is omitted. For simplicity, we chose the size parameter b for the (Os)^ dineutron cluster as 
b = l/\/2v, where v is the width parameter ;y(^He) and j/(^He) optimized for the ^He and ^He, respectively, in the 
AMD calculations. The values v, which are listed in table HI correspond to hj. — 2.0 — 2.3 fm and satisfy the typical 
nn distance of the dineutron correlation. 



B. dineutron-cluster motion 

In order to investigate features of dineutron cluster structure in the 0+ states of *He, we extracted the ^n-cluster 
motion from the obtained *IIe(0+) wave functions. We assume a simple core ^He -f ^ n)o+ which is equivalent to the 
SU(3)-limit ^lle(O^), and form the 6He'^^(3)(o+)_2„ cluster wave function with the L — Q relative motion between 
the core ^He'^'^*^'^^ (0+) and the cluster. In the same way as [sj, [5l| for a-cluster motion, we calculated the 
reduced width amplitudes ry{r) for the ^n-cluster motion and the cluster probability S^^'^ by taking the overlap of 
the 6He'^'^(3)(o+)_2„ duster wave functions with the ^He wave functions. In Fig. [H we show the reduced width 
amplitudes in the *IIe(0^) and the ^lle(O^) wave functions obtained by the v58 and m56 interactions. These indicate 
the ^He'^'^^'^) (0+)-^n relative motion. We also show the reduced width amplitudes for the ^Be'^'^'^^^ (0+) — a relative 
motion in the ^^C(O^) and ^^C(Oj) given in 42]. Surprisingly, the ^n-cluster motion in the ^He is quite similar to 
the a-cluster motion in the ^^C. 

First we discuss the features of the dineutron clustering in the state. The most striking thing is that the 
®He(Oj) state has the large amplitude of the dineutron cluster in the long distance region around r = 4 — 6 fm, which 
well corresponds to the peak position of the a-cluster motion in the ^^C(Oj). The enhancement of the ^n-cluster 
component at the long distance is more remarkable in the v58 results than the m56 results. The cluster probability of 
the ®He(Oj), which is defined by the integrated overlap with the ^He^^^'^^ (0+)-^n cluster wave functions, is S^^'^ = 0.50 
and 5'^'^'^ — 0.43 in the v58 and the m56 results. The larger development of the ^n clustering in the v58 results is 
considered to be because of the stronger n-n interaction in the v58 than the m56 interaction. It is very important that, 
even with the weaker n-n interactions of the m56, the ^n-cluster structure survives with the significant component 
in the *IIe(0^). Considering that the other cluster exists inside the ^He'^'^'^^^ (O"*") core, it is regarded that the 
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FIG. 4: Reduced width amplitudes ryi=o{r) for '^He^^f3)(0+)-^n in the *He(0+), and those for ^Re-^n in the ^He(0+). The 
*He(0+) and ''He(0+) wave functions are calculated by AMD-f-GCM with (a) the v58 and (b) the m56 interactions. The 
^He^*^^^' (0^) is written by a SU(3)-limit ''He-2n cluster state. The ''He cluster and the cluster are expressed by the (Os)* 
and (Os)^ wave functions, respectively, where the size parameter for the (Os) state is chosen to be the same value as the 
AMD+GCM wave functions; b = l/y/2u{^Re) in the calculation of ry{r) for '^He(0+) and b = l/,/2i/(8He) in the calculation 
of ry{r) for ®He(0'''). The reduced width amplitudes for '*He-^n in the ^He^'^''^'(0"^) are also shown, (c) The reduced width 
amplitudes for ^Be-^^^^^ (0+)-a in the ^^C(0+) taken from [sj. 



®He(Oj) has the component of the developed "^Ke+^n+^n clustering, where two dineutrons are moving in L = 
orbits. Furthermore, from the analogy of the ^n-cluster structure in the ®He(Oj) with the a-cluster structure in the 
^^C, the ®He(Oj) is considered to contain the dineutron gas-hke structure. 

Next, we discuss dineutron structure in the ground state of ^He. In the ^IIe(0]^), the reduced width amphtude has 
a peak at the distance less than 3 fm. It means that the spatial development of the cluster is not so remarkable 
as that of the ®IIe(0^). After discussing dineutron structure in the ^lle(O^), we shall compare it with the dineutron 
structure in the ^IIe(0]^). In Fig. 01 we show the reduced width amplitudes of the "'He-^rt cluster motion in the 
^He(O^) obtained by the present calculations, and that in the ^He"^'^*^'^) (0'^) given by the SU(3)-limit '*He-^rt state. 
Compared with the SU(3)-limit, the calculated ^He(0|') wave function has a long tail of dineutron structure at the 
surface. The ^n-cluster probability in the ^lle(O^) state is S^'^'^ — 0.91 and 0.84 in the v58 and the m56 calculations. 
This is consistent with the fraction, 0.92 and 0.87, of the S = component, which are estimated from (5^). The 
^n-cluster probability is reduced by the S' = 1 component because of the mixing of the (^3/2)^ state. The dineutron 
wave function in the inner region is similar to that of the SU(3)-limit ^He-^n state. In this region, we have better to 
call it the spin-zero 2n correlation(dineutron correlation) rather than the cluster, because the antisymmetrization 
effect is important there. 

Comparing the result of *IIe(0^ ) with that of ^IIe(0]^ ) , we found that the reduced width amplitude for the dineutron 
component is suppressed in the ^IIe(0]'"). This is because of the P3/2 sub-shell closure effect. As mentioned in the 
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previous section, the j-j coupling feature is more remarkable in the ®He(Of ) than the ^He(O^). However, the cluster 
probability of the ^He(Of ) is still significant as S^^"^ — 0.57 and 0.52 in the v58 and the m56 results, respectively. 
This probability dominantly originates in the SU(3)-limit '^He+^n+^n configuration, which is equivalent to the L-S 
coupling p-shell configuration. It means that the dincutron correlation is still important in the ®He(0^). This situation 
is quite similar to that of the ^^C(O^) which is the admixture of the closure and the SU(3)-limit 3a state. As a 
result of the L-S coupling feature due to the dincutron correlation, the ^He(0]'") state should contain the significant 
(^3/2)^(^1/2)^ contamination. This result is consistent with the experimental indication of the Pi/2 component in the 
^He ground state reported by the recent observations pTl [28| . As seen in Fig.[4l it is also interesting that the ^He(0]'') 
state has a tail of the ^n-cluster motion at the surface, though the tail is slight compared with the long tail in the 
^He(O^). In conclusion, the ^lle(O^) is the admixture of the closure and the L-S coupling p-shell configuration 
of neutrons with a small tail of the dineutron clustering. 

C. condensate wave function 

In the previous subsection, we discuss the ^n-cluster wave function by assuming the core (^He +^ 't-)o+ which is 
equivalent to the SU(3)-limit ^IIe(0+). In this description, one of the clusters is confined in the the core (''He+^n)o+ , 
and its relative wave function to the ^He is given by the Is orbit of the harmonic oscillator potential with the oscillator 
frequency lj = 8v/3. 

As shown in Fig. [H in this SU(3)-limit, the radial wave function of the ^n-cluster around the "'He remains in the 
inner region. In such the case, although the ^n-cluster is moving in the S wave, the ^n-cluster receives much effect 
of antisymmetrization from the ^He core and it does not necessarily indicate a gas-like state. In order to see more 
directly the ^n-cluster gas-like nature, where two ^n's are moving in S wave far from the the "^He core, we assumed 
the condensate wave function in the ''He +^ n +^ n system and calculated the overlap with the obtained ^He(0"'") 
wave functions. 

We define the condensate wave function by naturally extending the a condensate wave function proposed by 
Tohsaki et al.[l| as follows, 

vE'cond(i?) = "0 / n {'^'S, exp (^-i^i^l^^ I $Brink(Sc, Si, S2, • • • Sfc), (10) 

where rip is the normalization factor and $Brink(Sc, Si, S2, ■ • ■ S^) is the Brink wave function for the C -|-fc(^n)-cluster 
system consisting the core(C) and k dineutrons(^n) as, 

$Brink(Sc, Si, S2, • • • Sfc) = ^ {0C(Sc)/"(Si)/"(S2) • • • /"(Sfc)} . (11) 

Here, the wave function of the ith ^n, (j) "(Si)j is given by the (Os)^ state localized around Si. Sc is the mean position 
of the center of mass motion of the core, and is chosen to be Sc = — -|(Si -I- S2 + • • • + Sfc). In heavy limit of the core 
mass A, this wave function is equivalent to the dineutron condensate wave function proposed by Horiuchi,^. In the 
present calculation for ^He-|-^n-|-^n, the core C is ''He, and the number of clusters is k — 2. We assumed the (Os)^ 
state of the core wave function, (jj and adopted the common size parameter b = 1/ •y/(2i/(^He) for the ''He and 
clusters. In the practical calculations, the 6-dimensional integrals for the coordinates. Si and S2, are performed by 
taking mesh points on (^12, |Si|, IS2I) and the total-angular-momentum projection (612 is the angle between Si and 
S2). 

In Fig. [5l we show the squared overlap, |(^He|5'cond(S)) p, between the condensate wave function and the ^He 
wave functions obtained by AMD-t-GCM. The calculated values are plotted as a function B which indicates the size 
of the spatial distribution of clusters in the condensate wave function. The ^He(0]^) has the overlap, about 0.5, at 
B < 2 fm. The condensate wave function ^'cond(5) with such a small size B is almost equivalent to the SU(3)-limit 
''He-|-2n-f 2ri state. On the other hand, the ^He(O^) has the maximum overlap, about 0.5, at remarkably large size 
B = A — 5 fm. This is an strong indication of the dineutron gas- like component in the calculated ^He(O^). The 
dineutron gas-like feature is further enhanced in case of the v58 interaction than the m56 interaction. These results 
are consistent with the discussion of the ^n-cluster wave function in the previous subsection. 

V. SUMMARY 

We studied the structure of ^He with a method of AMD-I-GCM. We chose the effective nuclear interactions by 
taking care of energies of subsystems, and reproduced the properties of ground states of ^He, ^He and ^He. In the 
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FIG. 5: The squared overlap between the dineutron condensate wave function ^'cond(-B) and the obtained *He(0"'") wave 
functions. See details in the text. 

ground state of ^He, the component of the P3/2 sub-shell closure is dominant. However, the L-S coupling feature is 
also significantly contained because of the spin-zero dineutron correlation. This is consistent with the experimental 
report on the significant mixing of {pz/2)^{Pi/2f component in the *He(0^). It is concluded that the *He(0]'') is 
the admixture of P3/2 sub-shell closure and L-S coupling p-shell configurations with a slight dineutron tail at the 
surface. This result is also consistent with the experimentally suggested large spectroscopic factor of the ^He(2+) in 
the SHe(0^'). 

The present results suggest that the state may appear a few MeV above the 2^ state. By analyzing dineutron 

structure, it was found that this state has a significant component of the developed ^He+^n+^n structure where two 
dineutrons are moving around the ^He core in S wave with a dilute density. The ■^n-clustcr wave function of the 
^He(02') state is similar to the a-cluster wave function of the ^^C(Oj) state. Therefore, we consider that the predicted 
0^ state is the candidate of the dineutron gas-like state, which is analogy to the a condensate state suggested in the 
^^C(O^). In the experimental energy spectra of ^He, some excited states were observed above the 2f state. Spins and 
parities of these states have not been definitely assigned yet. Since the present calculations predicted the remarkable 
neutron matrix element for the monopole transitions ^He(O^) ^*He(02 ), we expect that the *He(Oj) might be excited 
in inelastic scattering on nuclear target. 

Since the AMD framework is regarded as a kind of bound state approximation because of the restricted model space, 
coupling with continuum states is not taken into account. In future study, widths of the excited states should be 
carefully investigated by taking into account the continuum coupling in order to confirm the stability of the resonances 
against particle decays. 

In the present work, the calculations were performed within the AMD model space by using effective interactions. 
We chose the interaction parameters by taking care of subsystem energies such as a-n, ^He as well as nucleon- 
nucleon systems. Although it is difficult to completely reproduce all of the subsystem energies with a unique effective 
interaction, we found the interaction which can reasonably reproduce the global feature of the subsystem energies. 
We here stress that the level structiire of the excited states is not sensitive to the adopted nuclear forces within the 
reasonable choice of effective interaction, though the excitation energy relative to the ground state depends on the 
interaction. It is also important that the dineutron structure of the ®He(0"'") states is qualitatively similar among 
four sets of interaction adopted in the present calculations. For further investigations of He isotopes, more extended 
calculations based on the realistic forces should be important as well as ah initio calculations. 
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